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The charm-quark mass is typically not so far from the cutoff aT x in lattice simulations. Its deter- 
minant may then potentially introduce large cutoff effects. We choose the 0(a)-improved Wilson 
formulation and compute the vacuum polarization effects in two rather different observables at 
one-loop order. One is the quark-antiquark static force and the other the Schrodinger functional 
coupling; in addition we investigate two more quantities resulting from the latter. In all the cases 
the lattice artifacts due to the charm-quark are small when compared to the gluonic effects. This 
indicates that the inclusion of charm-quarks in dynamical fermion simulations is typically not a 
problem. 
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1. Introduction 

Some lattice groups have already started investing effort in studies involving charm-quark 
vacuum polarization in Lattice QCD simulations [1], The main goal of these is to exclude notice- 
able corrections due to charm-quark loops. One, for instance, could think of the non-perturbative 
computation of the quark masses and the A-parameter in QCD. If in a such a computation we 
can only use the three-flavor theory then the connection to the four-flavor theory can be done just 
perturbatively. The ALPHA-collaboration has already started computations in Nf = 4 using as a 
first observable the running-coupling in the massless Schrodinger functional scheme [2]. However, 
in order to set the scale, observables in the massive theory have to be computed. A typical lattice 
spacing is a ~ (2GeV) _1 . Hence, the question of whether such computations are feasible is arising. 

Our main worry occurs from the fact that in this case the charm-quark mass in lattice units is 
as large as am c = 1/2. The main problem results from the findings of [3,4] according to which 
in the valence-quark sector of the 0(a)-improved theory, although for masses up to am c = 1/2 the 
cutoff effects are sizeable and obey approximately a quadratic behaviour in a, above this mass the 
Symanzik analysis and improvement of cutoff effects breaks down. 

The purpose of our work is to investigate how big the cutoff effects in charm-quark vacuum 
polarization effects are. We, therefore, expand perturbatively a few observables in the renormalized 
coupling and investigate the dependence of the first non-trivial expansion coefficient on the lattice 
spacing and the quark mass. We make use of our previous work with H. Panagopoulos [5] and that 
by Rainer Sommer and S. Sint [6] in order to extract the lattice cutoff effects due to fermions in the 
qq-force F(r) and the Schrodinger functional coupling respectively. Using the latter we defined 
three different observables, namely the step-scaling function £ and the renormalized quantities v 
and p; for these quantities we extract the cutoff effects. This script is a compact version of Ref. [7]. 

2. Lattice Formulation 

We begin by commenting briefly on the most important features of the 0(a)-improved the- 
ory. All the relevant information is described in detail in Ref. [8]. According to Wilson's lattice 
regularization the total action S = S g + 5/ is given by: 

W = prX>(p)tr{l-I/(p)} and SfprfM =a 4 ££ ^(D + m 0> /)^, (2.1) 

p X (=1 

where the gluonic part S g is expressed as a sum over all oriented plaquettes p with U(p) the path- 
ordered product of the gauge fields around p. We consider an infinite lattice and the weight factors 
w(p) are set to unity for now. The Dirac operator D reads: 

3 3 

^ = 5l{7^(V*+V M )-aV*V /i } + c sw f £ a, v F, v , (2.2) 

H=Q /i.v=0 

with V^ and V* the forward and backward covariant derivatives respectively and c sw = 1 + [g^j . 
One can perform a systematic investigation of lattice cutoff effects after renormalizing the theory; 
hence, the cutoff effects depend upon the renormalization conditions. We, therefore, first use a 
massless renormalization scheme with scale pi and then a massive scheme at scale pi in order to 
investigate cutoff effects in lattice perturbation theory. 

At one-loop order, the renormalized coupling and quark masses 1 are given by: 

1 At higher orders in perturbation theory, the sructure of the renormalization and the improvement is more compli- 
cated when all quark masses are considered. More information can be found in Ref. [9]. 
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g 2 (n)=§oZg{gl' a v) and (i»,-)r(m) = "V + O (g 2 ,) , (2.3) 

respectively in terms of the improved bare coupling and improved bare mass [6, 8] : 

£§ = gl 1 1 + 0.01200(2)g§a £ m 0J j + O (g 6 ) and m^ = m 0>i (l - \amo 4 ) . (2.4) 

For a complete 0(a) -improvement a la Symanzik one has to use the modified-bare coupling and 
modified-quark-mass in a massless renormalization scheme. 

3. A systematic way of investigating the cutoff effects 

It is usually normal to consider observables the perturbative expansion of which in the minimal- 
subtraction scheme gy^ is known in the continuum. We, therefore, renormalize in the MS scheme. 
We do so by first moving from the modified-bare coupling to the lattice-minimal-subtraction (lat) 
scheme and then to the MS scheme. The lattice-minimal-subtraction scheme is defined in such a 
way so that we subtract the logarithmic divergences in a order by order in perturbation theory. For 
the perturbative order we are interested in: 

«L(M) = <?o z iat (lo>«M) w i th z iat {go,an) = 1 - 2&o£§log(a/i) + O (fjj) , (3.1) 

where b = (llJV e /3 -2N//3) /{4-n) 2 . We use observables O that depend on a single length scale, 
let us say, r and on the masses (m,)R. These have the following perturbative expansion: 

O = g 2 + O^(z,a/r)g 4 + ..., (3.2) 

where z = (zi, ■ ■ ■ ,Zn ( ) and z% = («i)r • r. After moving from the bare to the modified-bare-coupling 
and then to the lattice-minimal-subtraction scheme, the expansion of the observable takes the form: 

= O coat (rn,z,gl t (n)) (\ + S {ril,z,gla{v),a/n) , (3-3) 

written as a continuum part and a lattice part 5o(rjU,z,g^ at (ju),a/r). The latter is expressed as: 

8 {rn,z,gl t (n),a/r) = 8$\rn,z,a/r) + B^ ) (rH,z,a/r)gl t (^) + .. . . (3.4) 

From the expression ( fO| ) 8o(rH,z,gh L (n),a/r = 0) = 0. We can now move to the MS by applying 
the finite scheme transformation [6, 10, 11]: 

lat.MS 

*Loo =4s(m)- V4w+°(4w)' ( 3 - 5 ) 

c iat,MS = __^ +2 .i35730074078457(2)A^ c -0.39574962(2)A^ f . (3.6) 

The expansion of the relative cutoff effects is written as: 



8 (rn,z,g^(n),a/r) = 



= 8 ( L, (rn,z,a/r)g 2 m (n) + .... (3.7) 

cont Sms>K)r 

Since the theory is renormalized by minimal-subtraction, these lattice artifacts are intrinsically per- 
turbative and have no non-perturbative extension. Nevertheless, one can use combinations of 8q 
from different observables in order to obtain the expansion coefficients of the non-perturbative cut- 
off effects. In the longer write-up [7] we demonstrate how by renormalizing in the Schrodinger 
functional coupling g| F (L) for massless quarks we obtain the non-perturbative lattice artifacts 



at small couplings. Their fermionic part is approximately the same as that of Eq. (3.7) i.e. ~ 
8^' f \rn,z,a/r). 

However, it is more convenient to introduce a massive renormalization scheme for non-perturbative 
computations in QCD with a charm-quark. This is due to the fact that the mass of the charm-quark 
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is larger than the typical QCD scale, thus, it has reduced vacuum polarization effects which is most 
efficiently implemented by using a finite mass renormalisation scheme. 

We can define a massive scheme with scale \i = 1/ro through a particular observable Oq\ 

g 2 m (H,m R ) = O = g 2 + O^(z,a/r)g A + ... . (3.8) 

It is straightforward to show [7] that for an observable O at a different length scale r the relative 
artifacts in the massive scheme are given as a combination of those in the massless: 

8oJr/ro,z,g 2 m {ll,mK),a/r) = 8^\rn,rm R ,a/r)-8^(rn,rom R ,a/ro) ^(jti,m R )+... . (3.9) 

4. The qq-static-Force. 

We start our investigation by considering the force F(r) = 4^V(r) between static quarks. If 
one uses on-axis potentials as we do here, the most natural choice to define the derivative would 
be: F na i ve (r na i ve ) = - [V(r,0,0) — V(r — a, 0,0)] with r na i ve = r— |. However, the force can also be 
defined in terms of an improved qq-separation rj according to which the force has no effects at 
tree-level order in perturbation theory [12], thus, F t ree('"/) = C F /4nr 2 with C F = (AT? - l)/2N c . The 
improved separations rj for on-axis potentials were calculated in [12]. It is worth mentioning that 
the non-perturbative force in the pure gauge theory when defined in terms of rj has much smaller 
lattice artifacts compared to r na i ve . The results presented below are expressed in terms of rj while 
we have in addition looked at all quantities for r na i ve with no any worth reporting changes; hence, 
from here on we have r = rj. 

The one-loop corrections to the force are calculated as described in [5]. The expression of 
the force in the MS scheme after setting the renormalization scale to the natural choice y, = 1 jr is 
given by: 

r C F a m (\/r) 



r 2 



(l+/ 1 (z,a/r)a M§ (l/r) + 0((4 § )} . (4.1) 

The term f\ is split into a gluonic (g) and fermionic (/) contributions such as: 

fx = /i>A) + £/i,/fe,a/V), (4-2) 

with continuum expressions [13] (/£ = 0.57721566 . . . ): 

/U(0) = v[-| + lfrs], (4-3) 



/i./feO) = i 



In 



/CO 



(4.4) 



The relative lattice artifacts of the force are defined as 



Ims,K)r 



4 1) (z,a/r)f| § (l/r) + ..., (4.5) 



wi 



th 2 : 



N f 

4xS£\z,a/r)=4xS£' g \a/r)+4n J £s£' f \z h a/r), (4.6) 

i=\ 

The gluonic contribution to the force at finite a/r is not our main concern. However, we need it 
for a complete picture on the size of the lattice cutoff effects. We, therefore, provide numbers for 
Sp (a/r) extracted from Ref. [14]. Concerning the fermionic piece, we perform a re-evaluation of 



9 (I) 9 

The pre-factor An on <5^ appears because we multiply with a = g /(4k). 
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Figure 1: The fermionic cutoff effects 4%8f' f \z,a/r) forz = (O).z = 1 (A)andz = 3(V) 



di (z,a/r) for the relevant bare-masses according to Ref. [5]. For more details on the procedure 
adopted to extract S P (z,a/r) see Ref. [7]. 

In order to obtain a first glimpse on the relevant size of the lattice artifacts let us view the glu- 
onic part An8 F hg) (a/r). For r/a = 2.277, 3.3 12 and 4.319, 4nd F l ' g) (a/r) = -0.232(6), -0. 190(19) 
and —0.151 (42) respectively. In Fig. [T] we display results for the fermionic contribution to the lat- 
tice artifacts 4n8 F (z,a/r). The scale of the y-axis is about a factor two smaller than the pure 
gauge artifacts. Hence, the absolute value of the fermionic piece of the cutoff effects for a single 
flavor is much less than the gluonic. Concerning the mass dependence, we see that the fermionic 
cutoff effects depend very little on the mass. 

5. Schrodinger Functional 

We call Schrodinger functional the field theory in a finite space-time with volume L 4 , Dirichlet 
boundary conditions along the time direction and periodic in space up to a phase 6 for the quark 
fields. According to Monte-Carlo simulations [6] the advantageous value of 6 is 6 = n/5 while 
6 = is the natural, and more aesthetically appealing, alternative. The action has the form ( 24 ) 



with all fields being zero outside < xq < L, the gauge fields having fixed values at the bound- 
aries xo = and Lq = and the fermionic fields being zero at these boundaries. For a complete 
0(a) -improvement one has to modify the weight factor of the timelike plaquettes attached to the 
boundaries w(p) = c t = 1 - [0.08900 + 0.0191417Y f ±0.00005]^ + . . .; for more details we refer 
toRefs. [6,8,15]. 

For our perturbative calculation we use the definition of the coupling appearing in Refs. [6, 15]. 
The Schrodinger functional coupling depends upon L playing the role of the inverse renormaliza- 
tion scale, the phase 6 and the dimensionless parameter v which appears on the fixed boundary 
gauge fields and is usually set to zero. In a non-perturbative computation of the running of the 
coupling the main object needed is the step-scaling-f unction defined as: 

L(u,z,a/L) = gl F (2L,2z)\ g 2 F{Lz)=uXmi)RL=Zi =u + Z l {z,a/L)u 2 + ... . (5.1) 

However, in order to have a more general picture we can also look at other quantities such as: 

a(u,z,a/L) = v(L, z) |^ F(L)2 ) = „ ](mi .) RL=Zi = vi (z, a/L) + O(k) , (5.2) 

with the quantity v defined explicitly by: 

- 2{T \, r = - 2 /. r -Vv(L,z), (5.3) 

g 2 v(L,v,z) £sf( l > z ) 
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Figure 2: Upper-Left: The cutoff effects An&£' 8) (a/L) (▼) and An8 y z J, (z,a/L) forz = 0(«),z = 1 (A) 
and z = 3 (■) extracted for 9 = n/5. Upper-Right: An8^ ,f) (z,a/L) for z = (•), z = 1 (A) and z = 3 (■) 
extracted for = 0. Lower-Left: The cutoff effects A^' g \a/L) (▼) and A^' f) (z,a/L) forz = 0(» (0))and 
z = 3 (■ (D)) extracted for 9 = n/5 {9 = 0). Lower-Right: The cutoff effects 4nA p Lf) (z,a/L) forz = 1 (• 
(O)), z = 2 ( A (A)) and z = 3 (■ (□)) extracted for 9 = n/5 (9 = 0). 



and 



p(u,z,a/L) 



g| F (L,z)-gj F (L,0) 



#|fM) 



gl v (L,0)=u,(mj) R L=Zi 



P[(z,a/L)u + 0(u 2 



(5.4) 



For the step-scaling function £, we consider the relative lattice artifacts: 

s i n\ Z(u,z,a/L) -L(u,x.,0) m. 

oz(u,z, a/L) = — ^ / — nN = 8^ ; (z, a/L)u + ... , 



£(w,z,0) 



L(u,z,a/L)-L(u,z,0) _ (i) ( 
z,0) 
while for both p and v the absolute artifacts: 

Av(u,z,a/L) = £l(u,z,a/L) — O(k,z,0) = A T , '(z, a/L) + 0(u) , 

A p (u,z,a/L) = p(u,z,a/L) — p(u, z,0) =A p (z,a/L)u + 0(u ). 

The leading perturbative terms 8^ (z,a/L), A\-, '(z,a/L) and A p '(z,a/L) can be decomposed into 
gluonic and fermionic pieces: 



(5.5) 

(5.6) 
(5.7) 



:(hg) 



8^(z,a/L) = 8^ 8) (z,a/L) + ^8^ n ( Zi ,a/L) 



(i,/)/ 



i-l 



N t 



4 l) (z,a/L) = A^\a/L) + ^A^( Zi ,a/L), 

Nt 

A p l \z,a/L) = ^A p lJ \ Zi ,a/L). 



(5.8) 

(5.9) 

(5.10) 
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In the upper-left and upper-right plots of Fig. ||we show the artifact 4k8^ '(z,a/L) for G = 
n/5 and 6=0 respectively. The cutoff effects for an individual 0(a) -improved fermion are much 
smaller than those for the gluonic part. They do appear to grow with z, however, not very much. 
Of course one should always bear in mind that for more than one fermion flavor the cutoff effects 
add up accordingly. 

If we now move to the cutoff effects of the v which appears in the lower-left plot of Fig. |2[ we 
observe that are bigger compared to the step-scaling function (given that the overall magnitude is 
Vi ss 0.1). In contrast to what we observed before, here, the cutoff effects decrease with the mass. 

Finally, in the lower-right plot of Fig. |we present the absolute lattice artifacts in p. These 
have to be compared to the continnum values of \%p (z, 0) which range from 0.095 (0.086) for z = 1 
up to 0. 188 (0.170) for z = 3 at d = n/5 (6 = 0). 

6. Conclusion 

We have investigated the cutoff effects of several observables in two rather different setups 
from including an 0(a)-improved Wilson charm-quark at one-loop order in lattice perturbation 
theory. We always restrict ourselves to am c < 1/2. We observed that these are comparable to the 
cutoff effects of the pure gluonic part; for instance see the lower-left plot of Fig. ^ Asa matter 
of fact, in some cases the fermionic pieces of the cutoff effects are even smaller than their pure 
gluonic ones. This is demonstrated clearly in the upper-left plot of Fig. ||. The reason of studying 
these was to check whether they become very large for very-massive quarks and more specifically 
for masses associated with the charm-quark. In contrast to our expectations of seeing large cutoff 
effects, we observe that they do not grow much with the mass and sometimes even decrease for 
larger masses. The overall conclusion is that the lattice artifacts remain small as we increase the 
mass and, thus, the inclusion of charm-quarks in dynamical fermion simulations does not seem to 
be a problem with the current available lattice spacings. However, non-perturbative investigations 
of lattice artifacts have to be carried out. 
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